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Abstract. We investigate the cosmological dynamics of a class of supersymmetric chameleon models 
coupled to cold dark matter fermions. Supergravity corrections ensure that these models are efficiently 
screened in all astrophysical objects of interest, however this does not preclude the enhancement of 
gravity on linear cosmological scales. We analyse the background cosmology and solve the modified 
equations for the growth of cold dark matter density perturbations in closed form. Using this, we 
go on to derive the modified linear power spectrum which is characterised by two scales, the horizon 
size at matter-radiation equality and at the redshift when the chameleon reaches the minimum of 
its effective potential. The model includes a cosmological constant in the form of a Fayet-Illiopoulos 
term, which emerges at late times due to the coupling of the chameleon to two charged scalars. We 
examine the conditions under which this leads to viable background cosmology and go on to analyse 
the deviations from the ACDM predictions in the linear regime. We find that for reasonable values in 
the model's parameter space there is generically a region where the model's cosmology is viable and 
current measurements can be reproduced. A small discrepancy of the matter power spectrum from 
its ACDM counterpart can be obtained in a smaller subset of the parameter space. 
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1 Introduction 

The recent discovery of the acceleration of the expansion of the universe [1, 2] has raised some 
perplexing conundrums in cosmology. If the acceleration is due to a cosmological constant then why 
is its value so small compared to the predictions coming from particle physics? If the cosmological 
constant is absent, which still requires some mechanism to screen the contributions to the vacuum 
energy coming from the electroweak phase transition, then what is source or dark energy driving the 
expansion? One of the simplest alternatives is quintessence, a scalar field decoupled from matter 
slowly rolling down its potential such that its effective pressure is negative. This scenario comes with 
its own draw-backs: one must fine-tune the model parameters so that the acceleration begins around 
the current epoch and the complete decoupling from matter is unnatural. One generically needs a 
light scalar with mass of order Hq to account for dark energy and so any coupling to matter results in a 
new, long ranged fifth-force which would violate solar-system tests of gravity. This issue is not unique 
to quintessence, any model of dark energy introduces at least one new degree of freedom, however 
general relativity (GR) is the unique (up to Planck suppressed higher order curvature invariants) 
theory of a massless spin-2 particle [3] and so these dark energy theories are, in some sense, equivalent 
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to modified theories of gravity. This has prompted a recent interest in the subject [4], however the 
problem of fifth-forces remains. 

One way of avoiding this problem is to note that all of our current tests of GR have been carried 
out in our local neighbourhood and so there is nothing precluding large fifth-forces that are active over 
long ranges provided that there is some sort of screening mechanism whereby the modifications are 
screened out locally. Such a mechanism was found in the form of the chameleon mechanism [5, 6] and 
since then similar mechanisms such as the symmetron mechanism [7] , the environmentally-dependent 
D amour- Poly akov effect [8] have emerged, as has a second, independent mechanism, the Vainshtein 
mechanism [9], which is present in Galileon models [10] and massive gravity [11]. In this work we will 
focus exclusively on the chameleon mechanism, which are a sub-class of scalar-tensor theories. 

Chameleon theories have been well studied (see [12, 13] for some reviews), both in the context 
of observational signatures and their cosmological behaviour, however there has been little progress 
towards any sort of UV completion. Ultimately, one would like to realise these models within fun- 
damental physics like string theory and a supersymmetric extension of these theories would be the 
first step towards this goal. In fact, [14]^ have used the KKLT mechanism [17] to find a chameleon 
coming from type IIB string theory compactifications at the cost of using the opposite sign in the 
gaugino condensation superpotential from the canonical one, which may act to decompactify the ex- 
tra dimensions rather than screen fifth- forces [18]. By examining supergravity breaking in a hidden 
sector, which can result in a direct coupling between the dark and observable sectors, the authors of 
[19, 20] have studied the chameleon mechanism which can emerge and have found a no-go theorem 
[21]: the scale of supersymmetry breaking is so large that it renders any effects from the matter 
coupling negligible. This effect may be avoided if one considers only global supersymmetry where the 
dark and observable sectors are secluded. The authors of [22] have done just this in the context of 
the supersymmetron^ . 

Scalar-tensor theories are IR modifications of GR and so such a bottom-up approach is sensible. 
In a recent letter [23] we have presented a general framework for embedding scalar-tensor screening 
mechanisms into global supersymmetry and have derived some model-independent features that arise 
when such theories are supersymmetrised. Within this framework, the chameleon couples only to cold 
dark matter (CDM) so that the observable and dark sectors only interact weakly via the breaking 
sector. This allows us to circumvent problems arising from hierarchies between the electroweak and 
dark energy scales. The mere presence of underlying broken supergravity is enough to ensure that 
supersymmetric symmetrons cannot exist in supersymmetry. Chameleons and dilatons survive. How- 
ever, with the exception of those arising from no-scale Kahler potentials, they are always so efficiently 
screened that no astrophysical signatures (such as the constraints of [24-29]) are present. Supersym- 
metry is always broken at finite densities, however, the seclusion of the dark and observable sectors 
ensures that the scale of this breaking is set by the ambient density and the model parameters and is 
generally well below the TeV scale associated with supersymmetric particle physics. 

In this work we will generalise the supersymmetron to an extended class of chameleon theories 
using this framework and investigate their cosmological dynamics. After introducing the general 
models in section 2 and the supersymmetric models in section 3.1 we will study their cosmology in 
section 4. These models have locally run-away scalar potentials which terminate at a supersymmetric 
minimum. Like all chameleons [30] , these models require a cosmological constant in order to account 
for the present day acceleration and energy density in dark energy with the exception of a very 
fine-tuned case. Going beyond the background level we proceed to study the linear perturbations 
of the cold dark matter. Unlike previous works which have imposed constraints by demanding that 
the modifications to the perturbation equations are negligible, we can solve the modified equation 
for the time-evolution of the density contrast in closed from using modified Bessel functions. This 
allows us to calculate the CDM power spectrum analytically and investigate the new features due to 
modified gravity. In particular, when the field converges to its supersymmetric minimum sometime 
between matter-radiation equality and the present epoch there are three separate regimes where the 
spectrum exhibits different scale dependencies rather than two. We first investigate the case where 

^see [15] for a generalisation of the model and [16] for an application to inflation. 
•^The supersymmetron is a chameleon and not a symmetron. 
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no cosmological constant is needed and find that large deviations from the GR power spectrum are 
predicted and so we abandon this scenario and calculate the power spectrum for the general case. 
By demanding that this does not differ too much from the GR prediction we find constraints on the 
model parameters. 

Next, we turn our attention to the need for a cosmological constant, which is difficult to include 
in supersymmetry. Supersymmetry is broken when the vacuum energy is positive and so we cannot 
add one at the level of the action and the contribution from supergravity breaking is of order M^^m^^^ 
(TO3/2 is the gravitino mass), which is far too high and must be fine-tuned away. In the same letter 
[23] , we addressed this problem by introducing a mechanism by which a coupling of the chameleon to 
two scalars charged under a local U(l) symmetry can act to drive the mass of one of the scalars to 
positive values at late-times, restoring the symmetry and leaving only a Fayet-Illiopoulos (PI) term, 
which acts as a cosmological constant. We assume that the old cosmological constant problem in the 
observable and hidden sectors is resolved and fix the value of this FI term so that this cosmological 
constant matches the observed present-day dark energy density. Unlike scalar vacuum expectation 
values (VEVs), FI terms run at most logarithmically and so do not suffer from matter-loop corrections. 
Therefore, whilst this choice is completely arbitrary within our globally supersymmetric framework, 
if one can find some natural reason for this small value in a more UV complete theory then it will 
remain at the same order of magnitude, even at low energy scales. 

In section 5 we apply this mechanism to our specific model. We indeed find a cosmological 
constant at late times, however at early times there are corrections to the effective potential which 
compete with the coupling to matter and act to negate the chameleon dynamics. We search the 
model parameter space for regions where the model parameters assume sensible values (to be made 
more precise below), the corrections are negligible before they vanish and the cosmological constant 
emerges before last scattering so that it leaves an imprint on the CMB. We find that such a region 
is ubiquitous and therefore these models are generically viable. We also find a region of parameter 
space where there are small deviations from the GR CDM power spectrum. We conclude in section 
6. 



2 Chameleon Gravity 

The chameleon screening mechanism may arise from the following action. 



d Xy/^ 



(2.1) 



which describes a scalar field coupled minimally to the observable matter but non-minimally to 
cold dark matter Xi via the Weyl rescaled metric (jf^^ — A^{(j3)g^^] ^(0) is known as the coupling 
function. In general, one may wish to couple the field to the visible matter as well. However, for the 
reasons explained in the introduction and re-iterated below, in this work we shall confine ourselves 
to a dark matter coupling only. As it stands, the action 2.1 describes the theory in what is known 
as the Einstein fvamc^ where the Ricci scalar is found using the Einstein fvamc metric g^i/ but the 
non- minimal coupling results in dark matter particles following geodesies of g^,^, the Jordan frame 
metric, so that observers in the Einstein frame infer the presence of a fifth-force 

F.^§^V^; m^M,,"^, (2.2) 
Mpi dip 

where ip is the canonically normalised field dip = k{4>) d(f). The coupling function is generally taken 
to be of the form A{(j)) = 1 -f 0{4)/M) -I- . . . with <C M so that perturbations in both frames do not 
differ too greatly. A second consequence of the non-minimal coupling is the emergence of an effective 
potential for <j). The equations of motion are 

U^^^-^Mt, (2.3) 
dip TOpi 
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where T is the trace of the energy-momentum tensor for dark matter (T^"^ = —2/y/—g6Sc/Sgf^^). In 
fact, T^'^ is not covariantly conserved in this frame since the dark matter fluid can exchange energy 
with the scalar; it is the Jordan frame energy-momentum tensor which is conserved V^T^"^ — 0. The 
trace of the energy-momentum tensor is equal to the density p for pressureless fluids, however since 
this is not conserved, it is more convenient to work with the non-relativistically conserved quantity pc 
defined by p = A{(f>)pc, which obeys the standard continuity equation. In the remainder of this work 
we shall treat pc as the conserved matter density, in which case the equation of motion 2.3 defines an 
effective potential 

KffM = ^M+Pc(AM"l). (2.4) 

The chameleon mechanism [5, 6] arises from potentials such as these when V{ip) takes on a run-away 
form and A((p) is monotonically increasing such that Vcs{(p) has a density-dependent minimum. When 
large, over-dense objects such as galaxies or stars are embedded into large, low-density backgrounds 
(for example, the cosmological vacuum) the field will try to minimise its effective potential, the 
minimum of which lies at different field values inside and outside the object. Chameleons have the 
property that the effective mass of small oscillations about the minimum 

mcs{ip) = V^^ + PcA^^ (2.5) 

is an increasing function of the ambient density. In over-densities of length scale R one typically has 
rriosR ^ 1 so that the force is very short ranged and is therefore negligible. If the object is large 
enough that the field can reach this minimum then the force is negligible and the object is screened, 
if not the mass of the field is a small perturbation about the background and the object is said to be 
unscreened. For a more thorough review of chameleon screening see [6, 12, 24, 27, 31, 32]. 



3 Supersymmetric Chameleon Gravity 

In a recent letter [23], we have introduced a general framework for embedding screened modified 
gravity into supersymmetry and have presented a class of supersymmetric chameleons with a locally 
run-away potential terminating in a supersymmetric minimum at large field values. We refer the 
reader there for the specific details; here we shall only specialise to the class of chameleon models we 
are concerned with. 

3.1 Supersymmetric Chameleons 

The Kahler potential for $ is non-canonical, which is a requirement for it to give rise to a run-away 
potential, whilst the dark matter fields have a canonical normalisation 

The self-interacting part of the superpotential is 

^=^f^V^f^V (3-2) 

where $ = (p+^/29ip+ . . . contains a scalar cj) whose modulus ultimately plays the role of the chameleon 
and $-1- ^ 4>± + ^/29ip± + . . . are chiral superfields containing dark matter fermions V'±- Splitting the 
chameleon field as (f>{x) = |0|e*^ and identifying cj) = \(j>\ from hereon, one can minimise the angular 
field (this is done explicitly in appendix B where a coupling to two U(l) charged scalars, which we 
will introduce later, is also examined) and define the new quantities 

A'^(^^Y'\,\ M"+4=(^^)''''Ao"+^ ^„.i„=(f)"Af, n = 2{a-P) (3.3) 
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to find the F-term potential 



Vf(0) 



dW 



When (j) <^ 0min this reduces to the Ratra-Peebles potential 



(3.4) 



(3.5) 



which has been well studied in the context of dark energy [33] (although one should be aware that 
we have not yet canonically normalised our field) . At larger field values the potential has a minimum 



at 



where l^(0n 



and dW/ d$ — 0. Supersymmetry is therefore broken whenever 



The coupling function is found by considering the part of the superpotential containing the 
interactions of $ and $± 

g 



1 



m A' 



<5-l 



Wint = m 

which gives a chameleon dependent mass to the dark matter fermions 



C D 



9$ I 



(3.6) 



(3.7) 



When the dark matter condenses to a finite density pc — m{'ip^ip^) this term provides a density- 
dependent contribution to the scalar potential resulting in the scalar-tensor effective potential Vcs = 
V + Pc{A — 1). With the above choice of superpotential, the coupling function is 



A(0) = 1 + —^0' = 1 
toAs 



(3.8) 



As it stands, the field is not canonically normalised since the kinetic term in the Lagrangian reads 



c 



kin 



The normalised field is 



(fi = Ai 



so that the coupling function 3.8 becomes 

A{ip) = l + x 

and the effective potential is 

Kff(</.) = A4 1 



Ai 



4^ 

mm 



mA; 



Xpc 



<5-l 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



which is shown in figure 1. We may then find the coupling f3{ip): 

-1 s , 



xSMpi 

f3(Pmin 



1 



^6^ Mpi ( ip 



mm \ y^min 



(3.13) 
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Figure 1. The effective potential. 



The effective potential 3.12 is the effective low-energy potential for a scalar-tensor theory de- 
scribed in section 2 with the scalar coupled to dark matter via the coupling function A{ip). Minimising 
the effective potential we have'^ 



n + S n + 25 



^ J \ ^ J Poo 



where 



5x 

At zero density the field sits at the supersymmetric minimum ip — Lp-min where its mass is 



Poo = Pc(l + ^co)' - inlMl,Hl{l + z^f = (3.15) 



- ^ = If + ^oo)^ ('^) Hi (3.16) 

2/3Vmin 2/32 VVmin/ 

When pc ^ Pea the field minimum is moved to smaller values by the matter coupling term and 
supersymmetry is broken. This is what was found in the general case studied in [23] . Supersymmetry 
is therefore broken locally in our model depending on the ambient density and Poo- The scale of 
this breaking is then set by the cold dark matter density and the model parameters, however this 
is generally far lower than the TeV scale associated with particle physics in the observable sector. 
This is one advantage of decoupling the dark and observable sectors, the dark sector does not suffer 
from an unnatural hierarchy of scales set by standard model particles decoupling in the visible sector. 
Away from the supersymmetric minimum, the field's mass is 



2(n + /3) fpn^inV^^ n + 2/3 /(^„,i„V"^+' , 2{S - (3) p, /(^,nin^' " 



n \ ip / n \ (p / n poo \ f 

(3.17) 

Clearly nias^'p) > TOqo when ip < (ySmin and so these models are indeed chameleons with a mass at the 
minimum of the effective potential which depends on the matter density. 

3.2 Supergravity Corrections 

In our previous work [23] we have found that once supergravity corrections are accounted for, super- 
symmetric models such as these screen so efficiently that no astrophysical objects can be unscreened. 



^The careful reader may notice that taking the hmit 6 = 1 gives a different equation from that found in [34], which 
contains a typographical error. 
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In particular, one can show that the self-screening parameter^ 



Xo - ...ry^ < f:^V < 10-"^ (3.18) 



where is the gravitino mass which can vary from ©(eV) in gauge mediated supersymmetry 

breaking to 0{TeV) in gravity mediated scenarios. When xo is less than the surface Newtonian 
potential, $n, the object will be screened. However, no object in the universe has (/)n < 10^"^^ thereby 
precluding any astrophysical fifth-forces. In this subsection we will briefly apply the reasoning that 
led to this result to our specific model in order to derive some constraints which will be important in 
the subsequent analysis. 

The correction to the potential coming from supergravity breaking is [23] 

which competes with the density dependent term in the effective potential (3.12). Since we focus on 
the branch of the potential where (j> < 0,nin, this term can always be neglected provided that it is far 
less than the density dependent term when </> has converged to its supersymmetric minimum. This 
requires that the supergravity corrections are negligible at densities around poo so that 

The denominator is proportional to the supergravity contribution to the vacuum energy, which is 
typically very large and is usually fine-tuned away whereas the numerator is proportional to the 
vacuum energy when the supersymmetric minimum is reached (see below), which we expect to be 
well below this. This condition tells us that the supersymmetric minimum must be well below the 
Planck scale and is simply the statement that the matter coupling and fifth-force is a low-energy, IR 
phenomena. It will be useful to express this condition in the alternative form 

« 3x170(1 + z^f { — )\ (3.21) 

from which it is immediately evident that (/?min ^ -^^pi even when the gravitino mass is as low as the 
gauge mediated supersymmetry breaking value of 1 eV. Using (3.15) equation (3.20) can be recast as 
a condition on the normalised field's mass at the supersymmetric minimum 



mlo » ^"^3/2. (3.22) 



The field's mass is at least as large as the gravitino mass. 



4 Cosmology 

In this section we will examine the cosmology of these models with the aim of accounting for dark 
energy. We will ultimately find that a cosmological constant is required in order to match both the 
present day equation of state w and the energy density in dark energy (this is a generic feature of 
chameleon cosmologies [30]). 

^See [24, 27, 32] for a more detailed discussion of the self-screening parameter. 
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4.1 Background Cosmology 

Solving (3.14) for the minimum in the hmit of both large and small dark matter density we have 




(4.1) 

Pc < Po 



We can now find the contribution to the vacuum energy density 



s 




x(5+n) l'poa\"+^ ^ 

Voff (¥>) « <( —f^P'^ [l^) > Pc > Poo (4.2) 

^ XPc Pc < Poo 

and the mass of the field using (3.17) 

I 9(^^,,) / r, \ " + « 

(4.3) 

Finally, one can find the equation of state for the field w^p. In uncoupled quintessence models this is 
simply Pip/ p^ = —1 when the field is at its minimum. However, in scalar-tensor theories the coupling 
of the field to matter results in a non-conservation of the density which forces us to define w via 
Pipl Pip = — 3-ff(l + w) [31] and so instead one has 

^^-(n^'^' • (4-4) 

I 0, Pc < Poo 

In order to match this with current observations we would like w — 1 and clearly this can be 
achieved by taking the limits pc 3> Poo, <5 S> «. and imposing the condition 

xb(\ + z^f « 3n. (4.5) 

This corresponds to the case where Zqo < and the supersymmetric minimum has not been reached 
by the current epoch. Both n and 8 appear as indices (or a combination of indices) in a superpotential 
and so we would expect them to be of similar order; taking (5 ^ n is then tantamount to neglecting 
many lower order operators in the superpotential, making the model appear somewhat contrived. 
When these conditions are not met, a cosmological constant is required in order to account for the 
present-day dark energy observations. It is important to emphasise that the need for a cosmological 
constant is not a new feature of supersymmetric models but rather a generic feature of scalar-tensor 
screened modified gravity [30]. What is new however is that it is not so simple to add one by hand 
within a supersymmetric framework. Global supersymmetry is broken if (P^) ^ and so the addition 
of a cosmological constant to the system is non-trivial. We will return to this issue in section 5. 

4.2 Linear Perturbations 

These models are always efficiently screened on small scales, however the screening mechanism is 
inherently non-linear and so one must still worry about fifth-force effects on the growth of structure 
on linear scales. It is entirely possible that the linear power spectrum is modified greatly in these 
models and it is important to impose constraints on the parameters such that the results of linear 
perturbation theory do not differ too much from GR. The growth of linear perturbations in screened 
modified gravity has been well studied [31, 32, 35] and the linear density contrast be = bpd Pc in the 
conformal Newtonian Gauge evolves on sub-horizon scales according to 

8, + 2i74 - \n,{a)H^ f 1 + '^^y\, I <5c - 0. (4.6) 
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Given (3.22), we have 07713/2 S> 2.5 x lO^^a Mpc ^ and so on the scales of interest we are always in 
the limit k <C rriipa. In this limit we can linearise (4.6) to find 



2H5r 



1 + 2k- 



'off'- 



0. 



(4.7) 



Equation (4.7) will be our starting point in what follows. It has solutions that can be written in 
terms of modified Bessel functions (whose degree depends on the model parameters) and so we have 
provided a short introduction to these functions, including how equations of the form (4.7) can be 
transformed into the modified Bessel equation, in appendix A. The quantity 5^ is gauge dependent and 
so is not a physical observable. Previous works have put constraints on their models by looking for the 
parameters where the final term in (4.7) is small and the GR result is recovered. In fact, it is the linear 
power spectrum that is the physical observable and so we shall constrain the model parameters by 
calculating the modified power spectrum and demanding that any deviation from ACDM is negligible. 
To this end, we shall make the approximation that modes do not evolve in the radiation era or outside 
the horizon during the matter era. Additionally, we shall treat the change from radiation to matter 
domination as a sharp transition and assume that the field settles into its supersymmetric minimum 
instantaneously at pc = Poo, thereby ignoring the short-lived transition period when equation (3.14) 
has no closed-form solution and any oscillations around the minimum (these are small on account 
of the large mass of the field [31]). Our powerspectrum therefore exhibits sharp discontinuities at 
the scales which enter the horizon at matter-radiation equality and (as we shall see momentarily) at 
Zjxj, which would be found to be smooth curves had we solved the full equations numerically. This 
calculation is aimed at discerning the new physics and features associated with our model and we do 
not compare the results with observation. These approximations are therefore acceptable. 

Physically, the last term in (4.6) corresponds to a scale dependent enhancement of Newton's 
constant: 

G..ik) 2Pi,)^^ ^^^^^ 



G 



= 1 



1 



On large scales the screening is effective and Gcff ~ G whilst on smaller scales the full enhancement, 
GcS — G{1 + 2/3((p)^) is felt. One would therefore expect that there is some wavenumber fc, which 
we will calculate presently, below which the modes feel no significant fifth-forces and the GR power 
spectrum is recovered. 

In what follows we shall consider two distinct cases. In the previous subsection we found that 
we can account for dark energy without a cosmological constant by imposing Zoo < 0, pc 3> Poo and 
S ^ n and so we shall first investigate this case. Finding that we cannot impose these conditions 
without large deviations from the GR power spectrum we will go on to examine the general case 
Zca ^ where it is assumed that the field reaches its supersymmetric minimum sometime around the 
current epoch i.e. Zoo 10 although our treatment will be valid for Zoo < -Zeq- 

4.2.1 The Cosmologically Viable Case 

This is the case where Zoo and 6 ^ n. Assuming a matter dominated era we can use equations 
(3.13) and (3.17) in (4.7) to find 



;i]°+9a;(l + zoo)'(fc<o)' 



Following appendix A, the growing mode solution is 



aktQ 



where 



■2 - ^ ' ^ + 30c ] and 



t 

to, 

81a;(l + Zoo)3 
16 



(4.9) 



(4.10) 



(4.11) 
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This should be compared with the GR prediction 



6cit) = CgkW"; where n = + + ^fio. (4.12) 

For small x, we have I,y{x) ~ x"^ [l + 0{x'^)] to leading order (see appendix A) and noting that 
= 6n + 1 we can see that these expressions agree for small k. Given the solution (4.10) we are 
now in a position to calculate the power spectrum. We start by noting that the time at which a given 
mode crosses the horizon {k = 2TTaH) is 

and assume that the modes are constant during the radiation era and outside the horizon in the 
matter era. In this case, the contrast during the radiation era (and outside the horizon in the matter 
era) is given by the primordial fluctuations from inflation, S^. Modes that enter the horizon during 
the matter era, that is modes with k < k^q = O.Olh Mpc^^, will begin evolving according to (4.10) 

and so we have the boundary condition Sdtu) = CyyQ,(}i)t-^^^ li,\GktQ(tYiltQ)^l'^\, which allows us to 
find CMG(fc) and hence the power spectrum 



^ (4.14) 

3^ k ^ 



Modified Bessel functions diverge from their leading order GR behaviour very rapidly and so if the 
modified power spectrum is not to deviate from the GR prediction too greatly the arguments of 
both functions in (4.14) must be small enough such that the leading order behaviour is a good 
approximation, at least over the entire range of k where linear theory is valid. Since to > this 
is equivalent to demanding akto <C 1. When this is the case one can verify using equation (A. 5) in 
appendix A that the leading order dependence is k^, the same as predicted by GR. Alternatively, 
according to (4.8) we should recover the GR prediction whenever k <C mes/l3{ip). Using equations 
(3.13) and (3.17) one finds an equivalent condition up to an order unity coefficient. Either way, one 
must impose that this is small when k — 0.1ft. MPc~^, corresponding to the scales at which the GR 
power spectrum becomes non-linear. One then has 

x{l + z^f <^ lQ-%-\ (4.15) 

Alternatively, one can define the above-mentioned scale, fc, below which no gravitational enhancement 
is felt: 



1 

P ~ — ^ ~ — — (4.16) 

This constraint is enough to show that this case is unnatural. If we demand that we have = — 1 
and Veg{(p) ~ HqM^^ then we must satisfy (4.5) simultaneously and so (4.15) imposes 6 ^ 3 x lO'^n. 
Since 6 appears as an index in the superpotential such a large value seems highly unnatural. With 
this in mind, we will abandon this limit and proceed to study the general case where we allow n and S 
to vary independently and ip to converge to its supersymmetric minimum at some point in the recent 
epoch. 

4.2.2 The General Case 

We start by noting that when Zoo > the final term in (4.7) will exhibit a different time dependence 
after the field has converged to its supersymmetric minimum, so we must keep track of modes that 
enter the horizon before and after this and match the time evolution appropriately. We therefore 



- 10 - 



begin by solving (4.7) for the case where z < z^o so that (f k. tp^i^ and m^p « moo- Assuming a 
matter dominated epoch and defining k^o — 2TTa(too)H(too) to be the mode which enters the horizon 
when z = Zqo (ioo = ^o(l + ^oo) ■^^^) we have 



(4.17) 



Following appendix A we can again write down the solution in terms of modified Bessel functions. 
Since the final term decreases with increasing t the growing mode is the modified Bessel function of 
the second kind: 

'to' 



Ckto 



where 



2 _ 



27x6 



(4.18) 



(4.19) 



Next, we must find the solution when k > koa- In this case we have (using equations (3.13), (4.1) and 
(4.3)) 



+ 



-r!^ + 9— - 1 + Zoo^fc'tg - 
3 n + Q \to 



(4.20) 



the solution of which is 



Scit)^C^l'-{k)I, 



aktQ 



to 



4(5 + n 



(4.21) 



with 



_ 81xS{d + n){l + Zoo)^ 
~ {AS + n)2 



and r = 



4(5 



3((5- 



(4.22) 



We can use this to calculate the power spectrum in the general case. Modes that enter the horizon 
during the radiation dominated era (i.e. k > fcoq) are constant until matter radiation equality when 



they start growing according to equation (4.21). In this case we have Cj^'°°{k > kcq)Ii,[cr{tcq/tkY] = 
^c(icq)- On the other hand, modes which enter during matter domination are subject to the condi- 
tion C^Q°° {k < kcq)Iv [(j{ti{/tkY\ = 5c{t'H)- Thus, modes that enter the horizon before Zoo evolve 
according to 



m = si 



Tktn 



(¥)' 



/„ \akt 
(Tktn 



to y 

^ *0 / 



(4.23) 



Near Zqo, modes inside the horizon (whether they entered during the radiation or matter era) thus 
evolve according to the general form Sc{t) = D{k)t^^^^I^[<7kto{t/toY] where the form of D{k) varies 
depending on when the mode entered the horizon as detailed above. When z = z^o the field converges 
to its supersymmetric minimum and the evolution now proceeds according to equation (4.18) and we 
must again match the two solutions ai z — z^o so that (5c (too) = D{k)tJ^^I^[akto{t^/toy']. Modes 
that enter the horizon later than this simply evolve according to (4.18), matching at the time when 
they enter the horizon, in which case we have (5c (in) — {k^-^^^^ Kui[Ckto{to/t}iY^^]. This leaves 

the power spectrum taking on different functional forms in three different ranges of k (this is to be 



- 11 - 



contrasted with the two predicted in GR): 



Pik) = (\Sl 



i t2 



i't\^kt o(t^y] Kiicktp] 

K^lCktp] 
KZlcktoi 



<k < k 



cq 



(4.24) 



k <^ 



We are now in a position to find constraints on the model parameters such that this power 
spectrum does not differ significantly from the GR prediction. We know from our analysis in the 
previous subsection that this requires taking the argument of all modified Bessel functions of the 
first kind to be small (<C 1), however the second kind functions require more thought. As detailed 
in appendix A, these grow with decreasing argument and diverge as it approaches zero and so one 
may be concerned that taking the argument to be small is not the correct limit. In fact, Ki^{y) ^ 
y-u Qf^y2-u;-^ gQ jj-j i-jjjg ijjjjji; Qjjg j-^g^y nBglcct the higher order terms. Once again we shall 
do this by imposing conditions such that the arguments of each modified Bessel function are small 
enough that one can neglect all but the leading order behaviour. One can indeed check by expanding 
the functions according to (A. 5) that this leading order behaviour recovers the GR result. This hence 
requires that the arguments of both the modified Bessel functions of the first and second kind are 
small, which, noting that to is the largest time in the problem, leads to two independent constraints: 



xS{S + n) 
{A5 + nf{l + z^) 



(4.25) 



coming from the largest argument of modified Bessel function of the first kind, ktQ(t^/tQY when 
k = O.lh MPc"^ and 

— ^ « 10-' h-' (4.26) 

coming from the largest argument of the modified Bessel function of the second kind, (kto, when 
k = O.lh MPc"^ We must also worry that the factor of CktQ{to/tii)^ in the denominator when 
k < koo can, in theory, be larger than C^^o when k = O.lh MPc~^. This only occurs when Zoo ^ 10^, 
however, for completeness, when this is the case we must impose the additional, weaker constraint 

x6 



1 



6(27r)2 



(4.27) 



5 Hybrid Dark Energy 

In a previous paper [23] we have laid out the general framework for dynamically generating a small 
cosmological constant at late times. We stress that we have nothing to say about the cosmological 
constant in the matter or any other sectors; we will assume that the old cosmological constant problem 
is solved. A mechanism like this is necessary, both because we have found previously that this model 
cannot account for dark energy without a cosmological constant and because it is not possible to 
include a bare cosmological constant at the level of the action in supersymmetry. 

5.1 A Late-Time Cosmological Constant 

The mechanism works by introducing two new scalars Il± ^ tt± + . . . with charges ±q under a local 
U(l) gauge symmetry. These have the canonical Kahler potential 

K ^ nle^''^n+ + nle-29^n_, (5.1) 

where X is the U(l) vector multiplet containing the gauge field and couple to the chameleon via the 
superpotcntial 

M^^ =.g'$n+n_. (5.2) 
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This construction gives rise to a new structure for the F-term potential as well as a D-term potential 
for the fields 7r±: 

V^^\{qnl-qnl-e)\ (5.3) 

where we have included a Fayet-IUiopoulos term which will later play the role of the cosmological 
constant. The new scalar potential is far more complicated with the addition of these new fields but 
when (7r_) = it reduces to our original effective potential for the chameleon (3.12) plus an effective 
potential for 7r+: 

F(^+) = ^(g7r2 -e')'+g'Vvr^; (vr-)=0, (5.4) 

where in this expression we have set 7r+ — |7r+| and will continue to do so from hereon. In appendix 
B we minimise the entire global F- and D-term potentials with respect to the angular fields coming 
from TT-t and show that (7r_) = is indeed a stable minimum of the system. 

The mass of the charged scalar 7r_|_ (or equivalently twice the coefficient of the quadratic term 
in (5.4)) is mj^ — g''^(j)^ — (Z^^^- At early times the chameleon is small (^ <j)rmn) and this mass is 
negative. The U{1) symmetry is therefore broken ((7r+) ^ 0). However, as the cosmological field 
evolves towards its minimum this mass slowly increases until it reaches zero, restoring the symmetry 
so that (7r+) = 0. We would therefore expect ttj^ — Q m the late-time universe leaving us with the FI 
term, which plays the role of a cosmological constant. Indeed, minimising (5.4) with respect to 7r+ 
one finds 

„2„2 



where 



& (5.6) 
V .9 

and = A is equivalent to the statement rrij^^ ~ 0. When (7r+) — equation (5.4) reduces to 
y(7r+) = ^^/2 and so we shall set ^ ^ 10^"^ eV in order to match the present-day energy density 
in dark energy. There is no natural choice for this parameter within our globally supersymmetric 
framework and so this value is completely arbitrary. It is worth noting however that FI terms are 
largely robust to quantum corrections; when supersymmetry is unbroken they do not run and when 
this is not the case they run logarithmically at most [36]. Therefore, if one could find a natural 
mechanism by which a small FI term is present in a more UV complete theory, for example a suitable 
combination of two or more large mass scales, then its value at lower energy scales will remain at the 
same magnitude"*; the same is not true of scalar VEVs, which receive large corrections from heavy 
particle loops. 

5.2 The Model Parameter Space 

Given the above mechanism, it is prudent to examine the model parameter space to determine the 
viable regions where a cosmological constant can appear and the model conforms with all experimental 
tests of GR. Since we have argued [23] that there are no astrophysical effects of the modification of 
GR, the only constraints so far come from the requirement that the linear power spectrum does not 
differ significantly from that predicted by GR. In what follows we will rule out the regions of parameter 
space where these are violated but will bear in mind that the parameters bordering these can lead 
to deviations from GR. If there were no other constraints then any set of parameters which satisfied 
these constraints would give a viable model, though one must be more careful than this. Firstly, when 
(tt^) ^ (i.e. at early times) there are corrections to the chameleon potential which can act to alter 
its cosmological dynamics. Secondly, we must ensure that the cosmological constant has the correct 
properties to reproduce current observations. We require the cosmological constant to appear before 
the present epoch and a necessary condition for this is 



> A or equivalently ( ^) " > [i^^ ' ^- (5.7) 



5 



Here we arc assuming that the cosmological constant problem in the hidden and observable sectors is resolved. 
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This is an additional constraint that must be imposed on the model parameters. Furthermore, if our 
model is to produce the correct imprint on the CMB then the cosmological constant must be generated 
before last scattering. We shall do this by imposing that the cosmological density pA (given in (C.4)) 
at which the U(l) symmetry is restored is greater than 1 eV*. 

5.2.1 Corrections to the Scalar Potential 

At late times (defined by the time at which = A) we have a FI cosmological constant, but at earlier 
times the non-zero VEV of 7r+ induces corrections to the effective potential for (f>: 

^corr - ^-^^^ " 1^0^ (5.8) 

q 2q^ 

These corrections compete with the density-dependent term coming from the dark matter coupling 
and therefore act to negate the chameleon mechanism. When they are important, they lead to a 
new, density- independent minimum and since the magnitude of the density dependent term decreases 
as the dark matter redshifts away it is possible to have a scenario where the field gets stuck at the 
new minimum and the cosmological constant is never generated. At first glance, the reader may be 
concerned that the correction proportional to —0^ results in a potential that is unbounded from below, 
but this form of the potential is deceptive. If one were to consider allowing the field to run away down 
this potential then at some point we would be in a situation where 4> > A and these corrections are 
no longer present; what looks like an unbounded potential is in fact a hill in the global potential. 

There are several possible scenarios involving these corrections, which either allow or preclude 
the generation of a cosmological constant depending on the model parameters. If the corrections 
are negligible compared to the density dependent term throughout the entire time that < A then 
they are never important to the model dynamics and vanish once </> > A. If, on the other hand, the 
corrections are important before they vanish then their dynamics must be included. However, if (jj 
can still pass A then a cosmological constant can still be generated since the corrections vanish after 
A is passed. If the only important correction is the quadratic one then a minimum always develops 
and therefore the cosmological constant will only be generated if the field value at this minimum 
is larger than A. If either the quartic correction or both corrections simultaneously are important 
then the potential may or may not develop a minimum. If no minimum develops then the field will 
eventually pass A since the potential takes on a (locally) run-away form. If a minimum does develop 
then we again require the field value at this minimum be larger than A in order to generate a viable 
cosmology. The exact details of how one can determine which scenario is applicable to a certain choice 
of parameters and whether or not the dynamics are affected to the extent that the model is not viable 
are rather lengthy and cumbersome and so they are given in full detail in appendix C. Below we shall 
only present the resulting parameter space once every possible scenario is taken into account. 

5.2.2 Low-Energy Parameters 

In order to classify the parameter space into viable regions we will need to derive certain conditions on 
combinations of the model parameters and so it is important to know which parameters are fixed in 
terms of certain combinations of the others. It will be sufficient to examine the position of the minima 
and the values of (j> relative to A and at no point will we need to use the dynamics of (p. For this 
reason, we will work exclusively with the field (j) and not its canonically normalised counterpart since 
this avoids unnecessary powers of /3. We have already seen in section 3.1 that three of the underlying 
parameters A^ (i = 0, 1,2) combine to form two derived parameters M and A. What are observable 
however are the low-energy parameters n, J, /3, cc, /z, Zoo, which are either combinations of M and A 
or indices that appear in the low-energy effective potentials (3.8) and (5.4); /i is a combination of the 
underlying parameter A3 and the dark matter mass m. It will prove useful to introduce the parameter 

G = g'/^q. (5.9) 

A static analysis therefore probes the six dimensional parameter space n,d,x, fi, Zac,,G and leaves (3 
unspecified. 
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In what follows, we will be interested in regions of parameter space where the cosmology is viable 
and the parameters themselves assume sensible values. In order to decide exactly what is meant by 
"sensible" it is instructive to pause and think about their physical significance, n and d are indices (or 
are combinations of indices) that appear in a superpotential and so these should naturally have values 
close to 1 as argued in section 4.1. g' = ^/qG is a U(l) coupling constant that appears in the coupling 
of the charged fields to and so we would expect values of 0(10^^ — 10^"^) so that the theory is not 
strongly coupled. Values much smaller than this would be tantamount to fine-tuning. Similarly, x 
parametrises the ratio of the vacuum energy density to the matter density when the field has converged 
to its supersymmetric minimum. The energy density due to the field today is (see equation (4.2)) 
^ff (0min) — xpc, which must be less than so that the dominant contribution to dark energy comes 
from the cosmological constant and so we require x < 0{1). Naively, one might argue that x should 
be small since it also parametrises the coupling to matter (see equations (3.11) and (3.13)) and so 
directly controls the enhanced gravitational force. However, we have already seen that supergravity 
corrections ensure that all astrophysical fifth-force effects are screened and so this argument does not 
apply. Linear perturbations are not sensitive to this efficient, non-linear screening and by demanding 
that their evolution does not differ too greatly from the GR predictions we obtained an upper bound 
on a; as a function of 2:00 in section 4.2. It is worth noting, however, that this is far weaker than the 
constraint we would have been forced to impose had the non-linear screening mechanism not been 
present ^. 

1 ^ 

Finally, we are left /i = m-sAg^ . When 6 = 1 this is simply the dark matter mass and thus 
varying /i is tantamount to fine-tuning the dark matter mass so that we get an acceptable cosmology. 
When (5^1 however we are free to fix the dark matter mass and what we are really varying is A3. 
In this sense we are not fine-tuning when we vary n but are in fact scanning the space of viable 
cosmologies as a function of A3. For this reason, we shall always fix 5 7^ 1. Now the dark matter mass 
can be any where from 0{eY) to ©(TeV) depending on the model and A3 appears as a mass scale 
in the underlying supersymmetric theory and so we would naturally expect it to be large (at least 
compared to the scales involved in the low-energy dynamics). Hence, in what follows we will treat 
anywhere in the region 0{eV) < /i < 0{Mp\) as sensible. 

5.3 Viable Cosmologies 

Given the above considerations and the procedure for dealing with corrections to the effective potential 
in appendix C we are now in a position to explore the parameter space. 

We have performed a thorough investigation into the exact effects of varying each of the six 
parameters on the cosmology and can find a large region where the parameters are indeed sensible 
and the cosmology is viable. It is difficult to gain any insight from the equations since they are 
all heavily interdependent in a complicated fashion and a large number of plots can be misleading 
since they can change very abruptly when a single parameter is varied by a small amount. For these 
reasons, here we shall simply describe the effect of varying some of the more constrained and less 
interesting parameters and present only a few two-dimensional cross-sections once these have been 
fixed at sensible values. 

Let us begin with the indices. nJ and 6 have very similar effects: if their value is increased whilst 
fixing the other five parameters then the region of parameter space where the corrections can be 
neglected will increase. Being indices, these should not stray too far from 0(1) and so their effects are 
far less pronounced than the other parameters, some of which may vary over many orders of magnitude. 
Hence, from hereon we will fix n = (5 = 2. Next, note that the constraints (4.25) and (4.26) coming 
from linear perturbation theory (which are schematically of the form {x/{l + ^oo)^)^^^ ^ 10^"" for 
y = 1 , 3 and w = 7,8 respectively) are really a constraints on how close we can push x to 1 given a 
value of Zoo- When Zqo = this essentially demands that x <C 10~®, which is equivalent to a small 

®Such a constraint is found in [34] where the non-Hnear spherical collapse of overdensities was used to find x <^ 10^'^''. 
In fact, this spherical collapse model applies to the unscreened situation (the supergravity corrections were not taken 
into account) and does not allow for the effects of screening once the overdensity is compact enough. It does, however, 
serve to highlight the sort of constraints that would be natural in the absence of screening. 

^Note that n = 2(/3 — a) and so really only even values of n are allowed. 
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coupling to matter required for small deviations from GR. Ideally, any theory of modified gravity 
should include a non-negligible coupling to matter, which is screened by other mechanisms than fine- 
tuning and so a more sensible model would have larger values than this. This can be achieved by 
raising z^o as is shown in figure 2 where we plot the Zoo-log(a;) plane with ^ = 10'^ TeV and G = 10~^ 
corresponding to what we have argued above are sensible values. For the sake of brevity we will set 
Zoo — 5 from hereon. This choice is completely arbitrary and different choices may give rise to very 
different cross-sections of parameter space, however, the region where the corrections are negligible is 
both ubiquitous and generically large when Zoo ^ and so it is not necessary to scan this parameter 
in great detail in order to narrow down a viable region. 




-20 -15 -10 -5 

Log(x) 



Figure 2. The various regions in the Zao-log{x) plane with n = S — 2, fi — 10'' TeV and G — 10~^. The yellow 
region shows the parameter range where the corrections are negligible. The magenta region shows the ranges 
where the quadratic correction is important, the dark blue region where both corrections are important and the 
grey region corresponds to the region where the constraints from linear perturbation theory are not satisfied. 
The zone at the boundary between the region where all constraints are satisfied and where corrections to the 
linear power spectrum are too large correspond to situations where a small deviation of the power spectrum 
from GR is present. 

Next, we plot the \og{fi)-\og{x) plane with n ~ 5 ^ 2 and G = 10^^ in figure 3 to investigate 
the effects of varying /i on the viable region. It is evident from the figure that large ( > ©(TeV)) 
values of /i are required for there to be a large region with negligible corrections; in fact, if one 
steadily increases /i one finds that this region grows, replacing the regions where the corrections are 
important. This behaviour is traced back to equations (C.4) and (C.7) in appendix C where it is shown 
that M^+" c>c /i" and therefore the density at which the corrections disappear increases slightly faster 
with /i than the densities at which the corrections become important. Finally, now that we have some 
idea of the viable values of Zqo and /i we plot the log(G')-log(a;) plane with n = S ^ 2 and /i = 10'^ TeV 



- 16 - 




-20 -15 -10 -5 

Log(x) 



Figure 3. The various regions in the log(/i)-log(a;) plane with n — S = 2, = 5 TeV and G — 10^^. The 
yellow region shows the parameter range where the corrections are negligible. The magenta region shows 
the ranges where the quadratic correction is important and the dark blue region where both corrections are 
important and the grey region corresponds to the region where the constraints from linear perturbation theory 
are not satisfied. The red region corresponds to models where a cosmological constant is generated after last 
scattering and are therefore excluded and the green region corresponds to models where ^min < A and a 
cosmological constant is only generated at some time in the future. The zone at the boundary between the 
region where all constraints are satisfied and where corrections to the linear power spectrum are too large 
correspond to situations where a small deviation of the power spectrum from GR is present. 

in figure 4 in order to investigate the values of G where the corrections are negligible. One can see 
that when G > C(l) the corrections are generally negligible, which is a result of (C.4) in appendix C, 
which show that the density at which the corrections disappear generally grows faster with G than the 
density at which they become important. The density at which the corrections are important both 
include an explicit factor of which is absent from the density at which the corrections disappear 
(there are other factors of x coming from the scale M though these vary with a far smaller power). 
This is the reason that the region where the corrections are negligible is larger when x is closer to 1. 
The plot clearly shows that there is a large region around G ~ 10~^, which we have argued above is a 
sensible range where there is no excessive fine-tuning or strong U(l) coupling. With the parameters 
we have chosen this only exists when x > 10^^°, however this does not really constrain x. If one 
were to increase either n, S, z^o or ^ this region would extend further in the direction of decreasing 
X. Thus, we have found that when the model parameters assume sensible values there is a large 
region of parameter space where the corrections to the effective potential are negligible before they 
vanish and an FI cosmological constant appears at late times. Clearly it is desireable to work in 
a regime where we do not need to concern ourselves with these corrections and so given that this 
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Figure 4. The various regions in the log(G)-log(3;) plane with n = 5 — 2, z^a ~ 5 TeV and /i = 10'' TeV. 
The yellow region shows the parameter range where the corrections are negligible. The magenta region shows 
the ranges where the quadratic correction is important and the dark blue region where both corrections are 
important and the grey region corresponds to the region where the constraints from linear perturbation theory 
are not satisfied. The red region corresponds to models where a cosmological constant is generated after last 
scattering and are therefore excluded and the green region corresponds to models where ^min < A and a 
cosmological constant is only generated at some time in the future. The zone at the boundary between the 
region where all constraints are satisfied and where corrections to the linear power spectrum are too large 
correspond to situations where a small deviation of the power spectrum from GR is present. 

region is ubiquitous at sensible parameter values we have not indicated in the sub-regions where the 
corrections are important but do not alter the dynamics. 

Before concluding, it is worth remarking that neither the linear perturbation analysis nor our 
analysis above has any constraining powers over the parameter /?. j3 is an index associated with the 
non-canonical kinetic term in the Lagrangian; if one works with the non-canonically normalised field 
(j) it appears in the equations of motion in the non-canonical factor ^(0) and if one instead works 
with if it appears in conjunction with n and 5 to give new indices in the effective potential. Since 
our analyses are performed using solutions of the equations of motion where Lp tracks the density- 
dependent minimum of the effective potential d^cj) ~ d^Lp ~ and it does not matter whether we 
consider a canonical kinetic term or not. This is the reason that /? is unconstrained. If one wishes 
to explore this parameter then it is necessary to go beyond this quasi-static approximation. Given 
this large field mass, this approximation is very accurate and so one would expect to find only weak 
constraints at best. 

Finally, let us notice that the only possibility of detecting deviations of these models from GR 
appears in the power spectrum and in the region of parameter space where all the cosmological 
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constraints are satisfied and the power spectrum is close, albeit different, from its GR counterpart. 
This occurs close to the boundary between the regions where the power spectrum differs drastically 
from GR and the one where no difference exists. As can be seen in figs 2-4, this is a non- negligible 
fraction of the whole parameter where new physics may be falsified. In the remaining parts of the 
parameter space, it seems that deviations from ACDM cosmology would be too minute to ever be 
noticed. 

6 Conclusions 

We have presented a class of globally supersymmetric chameleons coupled to dark matter fermions 
and have investigated their cosmological dynamics. Like all supersymmetric chameleons, except those 
arising from no-scale Kahler potentials, these models exhibit such efficient screening that there are 
no astrophysical signatures^. At small field values the scalar potential is locally run-away however, 
at larger values this behaviour terminates at a supersymmetric minimum. The coupling to fermions 
results in a density dependent effective potential and at finite CDM densities the field is displaced 
from this minimum to smaller values, breaking supersymmetry. Since the dark sector is secluded from 
the matter sector the scale of this breaking depends only on the model parameters and the ambient 
density so the scalar field's VEV does not receive TEV-scale corrections coming from the observable 
sector. 

The mass of the field is at least as large as the gravitino mass and so the density-dependent 
minimum of this effective potential is a stable attractor, which the cosmological field tracks from small 
values at early times to the supersymmetric minimum at late times. These models, like all chameleons, 
require a cosmological constant in order to match both the observed energy density in dark energy 
and the equation of state parameter. We have examined the growth of linear structures in these 
theories and, by solving the modified equation for the time-evolution of the density contrast in closed 
form, have derived the modified CDM power spectrum. Interestingly, this modified power spectrum 
exhibits three distinct scale variations rather than the two predicted by GR. By demanding that the 
modified spectrum does not differ significantly we have placed constraints on certain combinations of 
the model parameters. 

In order to address the issue of a cosmological constant, we have made use of a mechanism 
introduced in a previous letter [23] . A coupling of the chameleon to two scalars charged under a local 
U(l) symmetry results in a chameleon-dependent effective mass for one of the scalars. At early times 
this mass is negative, however as the cosmological chameleon evolves to larger values this steadily 
increases until it reaches zero, restoring the symmetry. When this happens, only a Fayet-Illiopoulos 
term remains in the scalar potential, which acts as a cosmological constant. If we assume that the 
old cosmological constant problem in the observable sector is solved and tune this FI term we can 
reproduce the small cosmological constant needed to match current observations. Our low-energy 
description does not include a natural value for this FI term, indeed it is arbitrary, however it is 
robust against quantum corrections and so if one could find a more UV complete model where there 
is a natural mechanism explaining this tuning then this value would survive to arbitrarily low energy 
scales. 

In the early universe, the symmetry is broken and this induces corrections to the effective po- 
tential. These act to negate the chameleon dynamics and so we have explored the model's parameter 
space in order to investigate whether this is an issue or not. We can find a large region where the 
model parameters assume sensible values, the constraints from linear theory are satisfied and the 
cosmological constant appears before last scattering. Thus we conclude that these models have a 
generically viable cosmology. 

Finally, whereas we can find a large number of viable models, there is no method by which they 
can be distinguished observationally for large parts of the parameter space. The efficient screening 
precludes all astrophysical signatures and the cosmology is indistinguishable from ACDM. We have 
seven free parameters, all of which can be tuned to sensible values to give a viable cosmology, one 
of which does not effect the linear perturbations or the corrections to the potential owing to the 

®We have not investigated a coupling to photons, which is not bound by these general results. 
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quasi-static motion of the field. On the other hand, we have also found a region in parameter space 
where small deviations of the power spectrum are present and therefore the models could be tested. 
It would be useful to investigate this further and look for mechanisms such as a coupling to photons 
or the MSSM which may help to break this degeneracy and provide more observational handles 
on these models. This is a daunting task since any coupling to the observable sector brings with 
it a plethora of complications. This is left for future work and here we shall be content that we 
can find supersymmetric models of screened modified gravity with viable cosmologies, despite their 
observational draw-backs. 
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A Modified Bessel Equations 

We have seen in section 4.2 that the linearised perturbation equations are all of the form 

t^Sc + atSc - + c^t^"-) Sc = 0, (A.l) 

with a = 4/3. The substitution 6c — f^Sdt) with n = (1 — a)/2 may be used to find the following 
equation for 6c'. 

t^!c + tic-( ^-^—^ +b^ + chA 6~c = 0, (A.2) 



(a -1)2 52 



4 

which may further be transformed into the form 



u^6c" + u6c' - {ly^ + u^) 6c; = ^-^^ + ^ (A.3) 

using the substitution u — df /r and the notation / = d/du. Equation (A. 3) is a modified Bessel 
equation^ the solutions of which are modified Bessel functions of the first and second kind, I^{u) and 
K^iu). Unlike regular Bessel, functions which are oscillatory in nature, these functions either grow 
{Ijj) or decay {Ki,) with increasing u. The general solution is then 



6c{t) = [Cih [-/) + C2K, (^f)] (A.4) 



although in section 4.2 we shall only be interested in whichever function is the growing mode (this 
depends on the sign of r). The modified Bessel function of the first kind has the power series expansion 

where r(TO) is the gamma function. The modified Bessel function of the second kind is defined via 

K^{x) ^ Imi — — — ^ . (A.6 

n^v 2sm(n7r) 

Its power series expansion is complicated and includes both negative and positive powers of however 
when X the leading order term is 

^"^("^^^ 2I^''"'+^('^'")- (^-^^ 



- 20 - 



B Minimisation of the Global Potential 



In this appendix we show how the global F-term scalar potential can be minimised to eliminate both 
the angular fields and 7r_ to recover the simple form given in equation (3.4). Ignoring the contribution 



from the Kahler potential for now (it depends on only) and setting 
minimum, we have: 



= 0, which is always a 



dW 





' |,/.|2"-2 


















At' 



/3 5ne 



ApAf^ 



(B.l) 



which, as can be seen, simplifies greatly when the negatively charged field has zero VEV. We shall 
see now that this VEV does indeed minimise the potential. We begin by writing the charged fields in 
polar form n± = Tr±e*^±. Ideally, one would hope to set the three angular fields {9,9±} to constant 
values in order to give negative signs in from of the final three terms in (B.l)^, however this is not 
possible and instead one must eliminate them in terms of the other fields. In order to do this, we 
exploit the local U(l) symmetry, which acts as 0±{x) — > ±qa{x), to set 9+{x) = 9^{x) = x(2;)/2, 
which reduces the angular fields to the set {9,x}- With this in mind, the scalar potential, including 
the contribution from the Kahler potential is 



2g'\Wi /Ai 



y(^,6',7r+,7r_,x) = 

o / A \ 2/3-2 

9' Ml 



+ 



/32 



2/3-2 



Ai 



2^-2 




+ 



V2I3 



2/3-2 



A a-3 A P 
^^0 ^^2 



7r-|_7r_ 



cos[(a - 13)9] 

33 C0S[X - I 



A a-3 
^^0 
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A: 



/3-3 



cos[x - (/3 - 1)9] 



(B.2) 



Minimising this with respect to x one finds 

sm[x - (a - 1) 



a: 



a-3 



13-3 



A: 



33 sin[x - (/3 - 1)9] = 0, 



(B.3) 



which may be used in the equation found by minimising equation (B.2) with respect to 9 to find a 
relation between sin[(Q! — /3)9] and sin[x — (/3 — 1)9)] (or equivalently sin[x — (a — 1)9]): 



sin[(a - m = |7r+^_ sin[x - (/? - 1)9]. 



(B.4) 



This may be used to eliminate x from the potential to find (with tp fixed at its minimum): 



T/(0,0,7r+,^_) 
'Ai 
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y2 A 2q — 6 
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(B.5) 



^This is the approach taken when the charged fields are absent [22] . 
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At first glance, one may worry about the square roots, however it is important to note that the 
above expression is only true when ijj is fixed to its minimising value. Furthermore, if one examines 
equation (B.4) then it is evident that as 7r-|_,7r_ — )• the second term in the square root has exactly 
the same behaviour and so there is never a region in configuration space where the argument is 
negative. Minimising (B.5) together with the potential coming from the D-term and the '!t± terms in 
the superpotential, 

^ I {qnl + -qnl + g'' <f {A + ^ (B.6) 

with respect to 7r_ one indeed finds that (7r_) = is a solution. If one expands the global potential 
around this minimum by setting 7r_ — > (tt-) + (57r_ then the coefficient of the Jtt^ term is 

q{e-q^l)+9'M''- (B.7) 

In theory, this can be negative, however we have not yet finished minimising the potential. In section 
5 we learnt that there are two possible solutions for 7r+ given by equation (5.5) when 7r_ = and so 
we should check that these are indeed stable minima of the global potential. The case where 7r+ = is 
clearly a minimum since the negative term vanishes. The second case gives the coefficient as 2(7'^|0p 
and so in either case the coefficient is positive and the stationary point is a stable minimum. When 
(7r_) = equation (B.4) gives sa\[{a— f3)9] = and hence cos[(a— /3)0] = —1. Making this substitution 
in (B.5) yields the far simpler form of the potential given in (3.4). 



Vd 







dW 


d7r+ 


+ 


d7r_ 



C D-term Corrections to the Effective Potential 

In this section we briefly show how one can explore the regions of the low-energy parameter space 
where the corrections to the effective potential coming from the U{1) symmetry breaking solution 
(equation (5.8)) may render the model not viable. In particular, we will explain how the sub-regions 
where ((> passes A and the cosmological constant is still generated can be discerned. The corrections 
to the effective potential are of the form 



Korr = - ^</'' (C.l) 



3 



and we shall make use of the definition 5.9 and for brevity define Z = {1 + Zqo) 
C.l Late Time Importance of the Corrections 

We can estimate the density at which each correction becomes important and we can no longer neglect 
them by equating each one with the magnitude of the density dependent term in turn. In this case, 
one finds that the field values (pi at which the order i corrections are important are 

^ ' and (C.2) 



xp 



mm 



lS-4 

min 



(C.3) 



respectively. Now we can always begin the cosmic evolution far enough in the past such that the 
density is large enough that the corrections are negligible, in which case the field evolves according to 
the background cosmology detailed in section 4.1. As the field evolves, the coefficient of the density 
dependent term becomes smaller and the corrections will eventually become important. If this occurs 
before the field rolls past A then we must correct the dynamics appropriately. If, on the other hand, 
this occurs after the field has passed A then these corrections will no longer be present and we can 
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neglect them completely. Using equations (4.1) and (3.10), we can estimate the densities pt at which 
4> = 4>i and the density pA at which ^ = A: 



The condition that the corrections can be neglected is then pA S> /C^. In this analysis we shall take 
"much greater than" to mean an order of magnitude i.e. pA ^ lOpi. 

C.1.1 Mass Scales 

One must be careful that the parameters above are inter-dependent and it is important to keep track 
of which are fixed by specific choices of others. In the analysis of sectin 5 we consider the low-energy 
model parameters {n, (5, x, to, G, Zqo} independent, which completely fixes the derived scales M and 
A via equation 3.3. We start by writing equations (3.11) and (3.15) in the form 

10-3eVy n ^ ' 

<^min = a;^M, (C.6) 
which can be combined using equation (3.3) to find 

n 

These relations can then be used to eliminate the quantities M and A in equation (C.4) in favour of 
the low-energy parameters. 

C.2 The Quadratic Correction 

When p2 ^ PA, Pi, Pc the eff'ective potential is 



Kff(</>)«A4 



This is minimised at field values satisfying 



1 - 



in 2 



. (C.8) 



n+2 / , \ ¥+2 ^9^2 ? 2 



^miii 



lA4 



(C.9) 



When GCa;i/V < we have (f> ~ 0min and so this case is still viable provided that (fimin > A. If the 
converse is true then the minimum lies at field values 



iA4 



(C.IO) 



and demanding that this is larger than A we find that the parameters must satisfy 

r+4G"a;?^" > nA^ (C.ll) 
in order for the cosmological constant to be generated. 
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C.3 The Quartic Correction 

When the quartic correction is important the potential takes the following form: 



1/(0) = AM 1 



The minimum, if it exists, is given by the solution of 



¥+4 



n+4 



n 



A 



(C.12) 



(C.13) 



and so the only possible solutions have > (/)min- This means that when this correction only is 
important the field will always pass A at some time and generate a cosmological canstant. With this 
in mind, one may wonder if the case 0min < A is allowed since in this case the field can still pass 
A if the minimum lies at large enough field values. This situation is highly unnatural since once 
the corrections vanish the field lies at values greater than t/fmin and will subsequently roll backwards, 
reintroducing the corrections. Hubble friction will eventually reduce the amplitude of the oscillations, 
however this leads to a situation that is highly fine-tuned and sensitive to the initial conditions and 
so we exclude it. 



C.4 Simultaneous Corrections 

When both corrections are simultaneously important the effective potential, including the matter 
coupling, is 

V{^)=A'(l-(^y] +.pf^y-^c^^ + G'e<P'- (C.14) 



When including the density term one should technically solve the entire dynamical system in terms 
of f and Pc{t), however we can glean all the information we need if we simply set pc = 1 eV^. As 
mentioned above, we require that the field rolls past A before Pc = 1 eV^ so that the cosmological 
constant is generated before last scattering and this will be the case if the minimum is located at 
field values greater than this by last scattering. Technically, this condition is not sufficient since it 
only guarantees that the minimum is located at values greater than A by last scattering, not that 
the field passes A by this time, however the large mass of the field ensures that this approximation is 
sensible. Far enough in the past the corrections are unimportant and field tracks its minimum owing 
to this large mass. Eventually we reach the epoch where all three terms become important and pass 
to the regime where the density dependent term is negligible. This transition is smooth and so given 
the large mass we expect that the field should simply remain fixed at the new, density-independent 
minimum and therefore the dynamics should not differ largely from the static analysis we will employ 
here. 

The minimum satisfies the equation 



mm 



n+2 



xSpc ( (f>- 



2-5 



2G^^^ - 2G 



4 J,2 



0, 



(C.15) 



which must be solved numerically for the minimum given a specific set of parameters and for the same 
reasons given in section C.3, we will impose 0niin > A. If this has no solutions then the potential is 
runaway near Pc — I eV* and the field will be able to pass A. When solutions exist the parameters 
where the minimum occurs at field values larger than A are viable and those where the converse is 
true are not. 
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